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^ ■ ABSTRACT 

^ : 

^\ , Aims. We present a high quality multiwavelength (from 0.3 to 8.0 /im) catalog of the large and deep area in the GOODS 
1 Southern Field covered by the deep near-IR observations obtained with the ESO VLT. 
Methods. The catalog is entirely based on public data: in our analysis, we have included the ^43514^, F&QQW , F775W and 
O F850LP ACS images, the JHKs VLT data, the Spitzer data provided by IRAC instrument (3.6, 4.5, 5.8 and 8.0 /.im), and 
• publicly available U-band data from the 2.2ESO and VLT-VIMOS. We describe in detail the procedures adopted to obtain 
this multiwavelength catalog. In particular, we developed a specific software for the accurate "PSF-matching" of space and 
> ground-based images of different resolution and depth (ConvPhot), of which we analyse performances and limitations. We have 
O l' included both z-selected, as well as isTs-selected objects, yielding a unique, self-consistent catalog. The largest fraction of the 
X , sample is 90% complete at z ~ 26 or Ks ~ 23.8 (AB scale). Finally, we cross-correlated our data with all the spectroscopic 
^ ■ catalogs available to date, assigning a spectroscopic redshift to more than 1000 sources. 

' Results. The final catalog is made up of 14847 objects, at least 72 of which are known stars, 68 are AGNs, and 928 galaxies 
, with spectroscopic redshift (668 galaxies with reliable redshift determination). We applied our photometric redshift code to 
IL" ' this data set, and the comparison with the spectroscopic sample shows that the quality of the resulting photometric redshifts is 
. 5^ ' excellent, with an average scatter of only 0.06. The full catalog, which we named GOODS-MUSIC (MUltiwavelength Southern 
, Infrared Catalog), including the spectroscopic information, is made publicly available, together with the software specifically 
5_j ■ designed to this end. 

a ' 

_ - I 
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1. Introduction 

The advent of the Hubble Space Telescope and of the mod- 
ern giant telescopes has opened a new era in observational 
cosmology, where galaxy evolution can be traced back to 
its early stages. In this context, deep multicolour imag- 
ing surveys are established as a powerful tool to access 
the population of faint galaxies with relatively high effi- 
ciency. An already classical technique l|Steidel et al. 19951 
IMadau et al. 1998jl is used to select high redshift galaxy 
candidates with sharp colour selection criteria, which has 
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been shown to be particularly efficient at z ~ 3 — 4 
IjSteidel et al. 1999|l . 

Other strategies have focused on surveys that sam- 
ple the whole spectral range from the U to the K band, 
enabling galaxy evolution to be followed on a wider 
range of redshifts, and that have triggered the devel- 
opment and use of photometric redshifts to extract in- 
formation about galaxies at any redshift. The Hubble 
Deep Field North (HDFN), and the subsequent Hubble 
Deep Field South (HDFS), are the most successful ex- 
amples of this approach, which has stimulated an im- 
pressively large number of scientific investigations (see 
[Ferguson , Dickinson fc Williams 20 00 for a review). 

The Great Observatories Origins Deep Survey 
(GOODS. IDickinson 200l1l represents the next generation 
of this kind of survey. It is based on observations of two 
separated fields of 10x15 arcmin each, centred on the 
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HDFN and Chandra Deep Field South (CDFS hereafter), 
respectively. Over a total area that is about 50 times the 
combined HDFN/S, it is the result of an impressive effort 
with the most advanced facilities. 

The available dataset includes ultradeep images from 
ACS on HST, from the mid-IR satellite Spitzer, from 
the ultraviolet satellite GALEX, from the X-ray satel- 
lites Chandra and XMM, as well as from a long num- 
ber of ground-based facilities (see IGiavalisco et al. 20041 
for a more detailed presentation). Spectroscopic follow- 
up is also available or ongoing from a list of sur- 
veys and facilities fC imatti et al. 20021 ICowie et al. 20041 
IVanzclla ct al. 20 05) . 

In particular, the GOODS southern pointing, located 
over the CDFS (centred on a = 03 : 32 : 30 5 = -27 : 48 : 
30, Gi acconi et al. 2000|l and encompassing the ACS Ultra 
Deep Field (UDF hereafter. iBeckwith et al. 2003|l and the 
GMASS^ ESO Large Program (Galaxy Mass Assembly ul- 
tradeep Spectroscopic Survey, PI Cimatti), has been the 
target of extensive observations with ESO telescopes, car- 
ried on in the spirit of public surveys. To date, a ma- 
jor effort has been devoted to deep imaging observations 
in the J, H, and Ks bands with the VLT-ISAAC in- 
frared imager. The final area covered by the Ks images 
approaches 140 arcmin^ to a la limiting depth of typi- 
cally 26-27 mag/arcsec^ (AB scale), making it a unique 
combination of depth and size in the near-IR. Additional 
observations in the U band have also been obtained, and 
a large spectroscopic follow-up is currently underway. 

This paper describes the multicolour catalog, which 
we named GOODS-MUSIC^ (MUltiwavelength Southern 
Infrared Catalog), that we obtained using the public imag- 
ing data available in the GOODS-CDFS region. It in- 
cludes two U images obtained with the ESO 2.2m tele- 
scope and one U band image from VLT-VIMOS, the ACS- 
HST images in four bands, the VLT-ISAAC J, H, and Ks 
bands as well as the Spitzer images in at 3.5, 4.5, 5.8, and 
8 /im. Most of these images have been made publicly avail- 
able in the coadded version by the GOODS team, while 
the U band data were retrieved in raw format and reduced 
by our team. 

From this dataset we obtained both a z-selected, as 
well as a iiTs-selected, sample. We also collected all the 
available spectroscopic information and cross-correlated 
the spectroscopic redshifts with our photometric catalog. 
For the unobserved fraction of the objects, we applied 
our photometric redshifts code to obtain well-calibrated 
photometric redshifts, which are described in a forth- 
coming paper (Vanzella et al. in preparation), together 
with the redshift based on the neural network technique 
IjVanzella et al. 2004|l . 

The resulting dataset will be used for several types of 
scientific output, for which we refer to forthcoming papers. 

1 ^http:// www. arcetri. astro, it/^ cimatti/gmass/gmass. html 
^ Note that our sample is definitely unrelated to 

the data and activities of the MUSYC collaboration, 

ihttp : / / www . astro . yale . edu/musyc/| 



At the same time, we make the dataset publicly available, 
along with some software specifically developed for this 
project, all of which can be of wider interest, and which 
we first describe here. 

The paper is organised as follows. In Sect. 2, we briefly 
summarise the available dataset. In Sect. 3, we discuss 
the procedure for object detection in the ACS-z image, 
and the corresponding measurement of total magnitude. 
In Sect. 4, we discuss the procedure adopted to estimate 
colours of z-selected objects, which has requiered the de- 
velopment of a new software code, which we explain in 
some detail. In Sect. 5, we focus on the i^s-selected sam- 
ple, describing how we identify the objects that are de- 
tected in Ks and not in z, and how we measure their 
total magnitude and colours. In Sect. 6 we briefly list the 
spectroscopic observations that we used to assign spectro- 
scopic redshifts, while Sect. 7 discusses the photometric 
redshift technique applied. Section 8 finally describes the 
format and the technical details adopted for the public 
release. 

All magnitudes are given in AB scale, and a A- 
dominated cosmology (fim = 0.3, f^A = 0.7 and Hq — 70 
km s^^Mpc^^) has been adopted throughout the paper. 

The data and photometric redshift cata- 
logs are also available on line at the WEB site 
Ihttp://! bc.oa-roma.inaiF.it/goods 

2. Imaging data 

2.1. GOODS-south ACS data 

In this paper we will focus on the southern field of the 
GOODS survey, located in the CDFS. 

Images are made available for all 4 GOODS filters, 
which are those used in the UDF observations: F435W, 
F606W, F775W, and F850LP. In the following, we will 
refer to these images as B, V, i, and z, respectively. Of 
this data set, in this work we use the data release 1.0 pro- 
vided by the GOODS team HGiavalisco et al. 2004|l . Table 
n summarises the characteristic of the each pointing of the 
GOODS ACS data. 

The VI. ACS GOODS images released by the 
GOODS team were divided into a grid of 4 by 5 8K x 
8K frames for each band. Since we were only interested in 
the portion of the image covered by the J and Ks obser- 
vations, we considered only the tiles that had an overlap 
(at least partial) with the IR images. Of these images, 
we built a single mosaic for each band (the software used 
for this purpose, named fitstile, is available on our web 
site): each resulting frame is 3.4Gb wide, which can still be 
handled by SExtractor in a single run, provided it is com- 
piled with the "large file" option. Although this solution 
is somewhat impractical (for instance, common visualisa- 
tion tools fail to display such large images), it allows us 
to avoid the complications arising when catalogs are ob- 
tained on overlapping images, such as inhomogeneities in 
the large scale background, different object detection, and 
other border effects. 
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2.2. The ISAAC data 

The GOODS field of the CDFS was the target of a 
deep imaging campaign in the near infrared with the 
ESO telescopes. A large field (20 by 20 arcmin) was 
covered with SOFI at a shallow magnitude limits in 
the J, H, and Ks bands. The GOODS-CDFS is being 
covered by much deeper observations in the same NIR 
bands with the ISAAC instrument. These data have 
been partially released by ESO and will be used in 
this work. A full documentation of this data set is pre- 
sented in Vandame et al. (in preparation), while details 
of the on-going GOODS program at ESO are given at 
\http ://www. eso.org/science/goods/products.html 

The version 1.0 data release, analysed with the 
ESO-MVM pipeline for imaging data reduction 
H Vandame 2002|l . includes 22 fully reduced VLT/ISAAC 
fields in J and Ks bands covering about 140 arcmin^ of 
the GOODS region in the CDFS. 

Figure ^ shows the position of the released ISAAC 
fields (yellow) superimposed on the optical image of the 
CDFS by ACS. The cyan square shows the position of 
the UDF on the GOODS-CDFS field, while the red rect- 
angle shows the layout of the K20 survey. The four white 
quadrants show the actual coverage of the VIMOS U band 
imaging. The H band observations cover only half of the 
entire survey, and are only limited to twelve fields. They 
have not been released by the present VI. data release 
but were part of a previous (V0.5) EIS data release and 
are public on the ESO archive. 

The typical exposure times range from 3 to 6 hours, 
the seeing ranges from 0.4 to 0.6 arcsec, and it is typically 
below 0.5 arcsec. We measured the value of the seeing for 
each ISAAC field by averaging all the bright stellar objects 
from the image, after normalisation of their profiles to unit 
fiux. For each ISAAC field, we also obtained a "convolu- 
tion kernel" /C, which is the transformation matrix that 
converts the ACS PSF to the relevant ISAAC one that we 
shall use for the colour estimate. This transformation was 
computed in Fourier space, as JC' = P^^isaacI ^^^ACS^ 
where the prime refers to the Fourier transform. An op- 
timal Wiener filter, which suppresses the high frequency 
fiuctuations, was applied in the Fourier domain to remove 
the effects of noise. 

2.3. U-band data: 2.2-WFI and VLT-VIMOS 

We searched the ESO public archives for deep imaging 
with the aim of complementing the spectral coverage by 
the GOODS team. Particularly important is the U band, 
which improves the photometric redshift estimate, in par- 
ticular for the lowest redshift {z < 0.5), when the U band 
still probes the blue side of the 4000 A break and, most 
important, is fundamental for identifying z ~ 3 galaxies 
(U-dropouts) . 

In the ESO Science Archive there are U-band im- 
ages taken with the wide field imager (WFI) at La 
Silla (Chile), which are part of the EIS public survey 




Fig. 1. The ACS GOODS-S field with superimposed the 
ISAAC tiling made by ESO to cover the whole field in J 
and Ks. The field shown is the ACS in the z-band limited 
to the area common to ACS and ISAAC observations. The 
cyan square marks the position of the ACS UDF. The 
NICMOS UDF Treasury observations cover a field that 
is inside the UDF. The red rectangle marks the position 
of the K20 survey. The four white quadrants show the 
actual coverage of VLT-VIMOS U band imaging; large 
gaps are visible, since the observing program has not yet 
been completed. 

IjArnouts et al. 2001|l . as well as recent images with the 
VLT-VIMOS imager. 

The WFI images were obtained in two filters, the so- 
called U35 and U3S. The Uss is the standard Bessel U with 
a peak efficiency of 50%, centred around Ac = 3800 A, 
while the U35 is a bit bluer, Ac = 3500 A, with higher 
efficiency (~ 80%), but with a red leak at A > 8000 A, 
where the WFI CCD is still sensitive. The net exposure 
time is 53654 s (- 15h) and 75100 s (- 21h) for the U35 
and C/38 filters, respectively. The U band image of VIMOS 
is based on a redder filter (Ac = 3900 A) and has an ex- 
posure time of 10000 seconds in total. The coverage of the 
GOODS-CDFS field, however, is partial, since the observ- 
ing program has not yet been completed. 

We decided to use these images even if they are 

affected by red leakage: by inserting the correct (with red 
leak) filter transmission curve to the photometric redshift 
procedure one is able to mimic the behaviour of distant 
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galaxies in the filter. The VIMOS images available on 
Dec 2004 cover a significant fraction of the field (about 
60%), although they exhibit the large gaps arising from 
the array disposition (see Fig. 

The raw U-band images were collected from the 
ESO Science Archive and reduced using Python scripts 
and the IRAF task mscred following the instruction of 



Valdes (2002) The total exposure times of these images 
surpass the previous data release of EIS-DPS, which only 
reached 43200 and 61200 seconds in the C/35 and Uss, re- 
spectively. The image quality of these images is good, with 
a seeing of 0.9-1.1 arcsec and a magnitude limit of 26.5- 
27.0 (AB mag) at 5 sigma in a photometric aperture of 
3 arcsec. The U image of VIMOS had good seeing con- 
ditions (0.7-0.8 arcsec) and reached a deeper magnitude 
limit (28 AB, 5cr). 



2.4. IRAC data from the Spitzer Legacy Program 

The GOODS survey also incorporates a Spitzer Space 
Telescope Legacy Program to carry out the deepest ob- 
servations with this facility at 3.6 to 24 microns, to study 
galaxy formation and evolution over a wide range of red- 
shift and cosmic lookback time. 

The first and second Spitzer data releases (DRl and 
DR2) consist of "best-effort" reductions of data taken with 
the infrared array camera (IRAC, IFazio et al. 2004|l on- 
board Spitzer. These are images from the two epochs of 
the "superdeep" IRAC observations for each of the two 
GOODS fields (Dickinson et al., in preparation). 

These fields were imaged at 3.6 — 8/i with IRAC, with 
a mean exposure time per position of approximately 23 
hours per band (doubled in the overlap strip containing 
the UDF), reaching far deeper flux limits than observa- 
tions planned for the Guaranteed Time programs. 

The IRAC observes simultaneously in all four chan- 
nels, with channels 1 and 3 (3.6 and 5.8 microns) covering 
one pointing on the sky, and channels 2 and 4 (4.5 and 8.0 
microns) covering another. The two IRAC flelds of view 
are separated by about 6.7 arcmin in the focal plane, and 
the long axes of the GOODS fields are oriented along the 
direction separating the two IRAC fields of view. The con- 
sequence of the 2x2 mapping pattern is that, in a given 
observing epoch, the area covered by channels 1-1-3 and the 
one covered by channels 2+4 have a small region of over- 
lap (about 3 arcmin). The GOODS-S four-channel overlap 
region includes the Hubble UDF IjBeckwith et al. 2003|l . 

The first and second epochs of IRAC images publicly 
released by the GOODS Team, cover the entire GOODS- 
CDFS field at 3.6, 4.5, 5.8 and 8.0 (see Fig.EJ. Tabled 
gives detailed information for the GOODS-CDFS survey 
used in this paper, with the wavelengths, area covered, 
and magnitude limits for all the filters available at this 
moment. 
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Fig. 2. The four channels of the IRAC images for the 
GOODS-South field. The pointings are chosen to have 
the UDF in the overlapping regions, where the exposure 
time is twice that in the external parts. The K20 layout 
is shown, as well. 



2.5. Accurate estimation of photometric errors 

The correct evaluation of the photometric errors for as- 
tronomical sources is bound to the possibility of having a 
realistic variance image associated with the measure im- 
age. It has become a common feature of the data reduction 
pipelines to produce "weight map" images that are usu- 
ally proportional to the exposure time of each pixel and, 
therefore, provide a measure of the relative S/N across the 
image. These images can be used by SExtractor to obtain 
a first map of the background r.m.s., which is the main 
source of noise for faint sources. However, the interpola- 
tions introduced by drizzling the images (shifting, rotat- 
ing, correcting distortion, and subsampling pixels onto a 
finer grid) result in correlations between pixels in the driz- 
zled science images. Therefore, the apparent r.m.s. back- 
ground noise that one measures in the image is smaller 
than the actual r.m.s., due to the effects of these cor- 
relations (for a more detailed discussion of weight map 
conventions and noise correlation in drizzling, please see 
Case rtano et al. 20011 especially Sect. 3.5 and Appendix 
A). 

To overcome this problem, we used a simple and effi- 
cient method to derive the r.m.s. from the scientific image 
itself, based on the correlation matrix of the pixels of the 
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Table 1. A summary of the photometric data of the GOODS-South field used in this work 



FILTER 




AA 


EXPTIME 


FWHM" 


PIXSCALE 


ZP 


AREA 


MAGLIM" 




A 


A 


s 


arcsec 


arcsec / px 


AB 


arcmin^ 


90% 


C/35 


3590 


222 


53654 


0.90 


0.23 


28.520 


143.2 


25.5 




3680 


170 


75100 


1.10 


0.23 


28.755 


143.2 


24.5 


UviMOS 


3780 


197 


10000 


0.80 


0.20 


32.500 


90.2 


26.5 


B (F435W) 


4330 


508 


7200 


0.12 


0.03 


25.65288 


143.2 


27.5 


V (F606W) 


5940 


1168 


6000 


0.12 


0.03 


26.49341 


143.2 


27.5 


i (F775W) 


7710 


710 


6000 


0.12 


0.03 


25.64053 


143.2 


26.5 


z (F850LP) 


8860 


554 


12000 


0.12 


0.03 


24.84315 


143.2 


26.0 


JlSAAC 


12550 


1499 


12600^^ 


0.45 


0.15 


26.000 


143.2 


24.5 


HiSAAC 


16560 


1479 


18000^= 


0.45" 


0.15 


26.000 


78.0 


24.3 


KsiSAAC 


21630 


1383 


23400^= 


0.45^= 


0.15 


26.000 


143.2 


23.8 


CHllRAC 


35620 


3797 


82800 


1.60 


0.60 


22.416 


143.2 


24.0 


CH2iRAc 


45120 


5043 


82800 


1.70 


0.60 


22.195 


143.2 


23.4 


CHSjRAC 


56860 


6846 


82800 


1.90 


0.60 


20.603 


143.2 


22.0 


CHAiRAc 


79360 


14797 


82800 


2.00 


0.60 


21.781 


143.2 


22.0 



a) For ground-based images, FWHM corresponds to the seeing, while for space-based data it is the PSF of the instrument. 

b) The limiting magnitudes (at 90 % completeness) are the mean values on the field, averaged on the positions and areas of 
every object. 

c) For J, H, and Ks filters, we report only the mean value of seeing and exposure time: detailed values can be found in 
Vandame et al. (in preparation). 



background. Indeed, as we explain better in Appendix A, 
the true r.m.s. of an image is given by: 

+m ^ N 
j— — m i—1 

where N is the number of pixels used for the r.m.s. calcula- 
tion. Hi the value of the i-th pixel on the source-subtracted 
image, and m the correlation length of the noise, typically 
of 10-20 pixels. 

In practice, we used SExtractor to convert the weight 
maps (available for each image of our data set) into the 
corresponding image of apparent r.m.s. We then estimated 
the correlation matrix in a homogeneous section of the im- 
age to compute the accurate r.m.s., and finally normalised 
the apparent r.m.s. image accordingly. 

We verified this procedure in the ACS images, where 
the weight images l|Giavalisco et al. 2004|l coincide with 
the expected inverse variance maps (i.e., l/c^) per pixel. 
In this case, we obtained an absolute r.m.s. of 0.001396, 
comparable to the value 0.0014 obtained with a traditional 
approach'^, as well as with the publicly available weight 
maps (i.e., l/cr^). 

For the U, ISAAC, and Spitzer images, we applied the 
same technique and found that the true r.m.s. is typically 
a factor 1.4 times higher than the r.m.s. that was directly 
measured on the science frames. As explained above, in 
the case of the ACS images, we used the available weight 
maps to compute the r.m.s. image by simply inverting the 
square root of the weight map. 

^ We polluted the z band images of GOODS with synthetic 
stars of different magnitudes and derived the true magnitude 
errors by the variance over 100 realisations 



2.6. Magnitude limits and effective areas 

The total area of the survey is the intersection of the im- 
ages in B, V, i, z, J, and Ks bands. The t/35 and U^g, images 
coves a much larger field of view (FoV), approximately 30 
by 30 arcmin, the IRAC images completely cover the field 
layout, while the VIMOS-U and the H band cover about 
half of the total field. Images in V, i, and z cover the 
same area, while B is smaller due to the lower number 
of dithering steps. Since we are interested in z and Ks- 
selected samples with good photometric redshift accuracy, 
we limit the sample to the area given by the intersection of 
the B and Ks images. In this sample, all the objects have 
UssUssBVizJKs and IRAC coverage, and about a half 
also have VIMOS-U and H band observations. The total 
area is 143.2 sq. arcmin. and the layout of the resulting 
field is shown in Fig. ^ 

Associating a limiting magnitude to this area is not 
trivial, since the exposure maps in optical and NIR are 
different and very complicated. We first decided to esti- 
mate an accurate map of the magnitude limits as a func- 
tion of position using r.m.s. images in z and Ks. The lim- 
iting magnitudes are computed at 1 sigma level and in 
an area of 1 sq. arcsec. Figures 13 and 01 show this in the 
GOODS well as the relevant distribution of the 

magnitude limit. Dark areas represent deeper exposures 
that correspond to the faintest magnitude limits. In the z 
band image we identified 6 main areas with a different in- 
terval in the magnitude limits; in Ks band the magnitude 
limit intervals identified are 6 (see Table EJ. 

This information was used as follows. We associated 
the 1 — a magnitude limit to each detected object (either 
in the z or in the if s-selected sample) as obtained from 
the "magnitude limit maps" at the corresponding posi- 
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Table 2. Magnitude intervals that we used to divide the 
overall sample into subsamples with well-defined mag- 
nitude limits and corresponding areas for the z and 
Ks bands. See also Figs. El and H MAGLIMl and 
MAGLIM2 (in magnitudes per sq. arcsec) define the lim- 
its of each magnitude bin. MAGTOT is the corresponding 
limit for the total magnitude of the objects, obtained using 
the simulations described in Sect. 3.2. 



FILTER MAGLIMl MAGLIM2 AREA MAGTOT 



1(7 la arcmin^ 

z 26.65 27.65 7.16667 24.65 

z 27.65 27.85 6.40742 25.65 

z 27.85 27.95 15.4276 25.85 

z 27.95 28.10 89.9600 25.95 

z 28.10 28.18 10.3641 26.10 

z 28.18 28.90 13.8741 26.18 

Ks 24.70 26.08 10.7698 21.60 

Ks 26.08 26.20 8.95936 22.98 

Ks 26.20 26.50 40.5733 23.10 

Ks 26.50 26.70 38.2161 23.40 

Ks 26.70 26.90 34.8572 23.60 

Ks 26.90 27.10 9.82067 23.80 



Table 3. SExtractor parameters for detection in the z 
band 



PARAMETER 


VALUE 


DETECT_MINAREA 


13 


DETECT.THRESH 


0.9707254 


ANALYSIS.THRESH 


0.9707254 


DEBLEND_NTHRESH 


32 


DEBLEND^INCONT 


0.05 


BACK_SIZE 


120 


BACK_FILTERSIZE 


9 


■ BACKPHOTO.THICK 


100 


CLEAN_PARAM 


3 



using the z band big mosaic as detection image. We used 
■an external flag image to eliminate the borders of the 
GOODS field. 

Because of the large number of tunable parameters, 
the use of SExtractor on a given data set is not straight- 
forward. After the first checks on the catalog, we realized 
that it is impossible to find a unique set of parameters 
to obtain an optimal detection over large FoVs, since in 
large and deep areas the objects have a wide range of di- 
mensions, going from the small and faint objects to the 
large, low-surface brightness but with high total luminos- 
ity near galaxies. In particular, the most critical parame- 
ters are those regulating the deblending and extension of 
the object (DETECT-MINAREA, DETECT-THRESH, 
DEBLEND-MINCONT, and DEBLEND-NTHRESH). 

To provide an optimal detection technique through- 
out all the GOODS field, we modified the SExtractor 
code, following our experience on the HDFS, described by 
Vanzella et al. (2001). The idea is to adopt a general set of 
SExtractor parameters for the global image and to define 
small areas (subframes) where a set of ad hoc parameters 
are adopted to optimise the detection of problematic ob- 
jects. In this case we adopted a set of parameters tuned 
to detect the faint, compact objects to obtain the bulk of 
the catalog. With these parameters, SExtractor typically 
oversplits bright, extended, and irregular sources (such as 
face-on spirals or low-surface brightness objects) in many 
fragments. To detect these objects, we obtained another 
version of the catalog using an extreme set of parameters 
(i.e. requesting very large area and no deblending) and 
visually inspected all the objects with magnitude z < 23 
and all the objects with spectroscopic redshift. In about 
15% of the inspected objects we decided to adopt tailored 
SExtractor parameters in a corresponding "subframe" . An 
investigation at fainter magnitudes has shown that such 
cases become very rare at z > 23, so we ignored them. 

We also slightly modified the SExtractor code to en- 
sure that the "back-annulus" is at least 1 arcsec wide (as 
suggested by the GOODS team) and to correct a bug for 
the estimation of the isophotal-corrected magnitude in the 
"dual image" mode when the r.m.s. of detection and mea- 
sure image differ by several orders of magnitude (see also 
the SExtractor mailing list). 



tion. The simulations described in Sect 3. made it possi- 
ble to translate this 1 — a magnitude limit into a com- 
pleteness limit for the object detection. Each catalog was 
therefore separated according to the magnitude limit, us- 
ing the bins listed in Table |3 As a result, each catalog is 
actually the combination of 6 (for z and for Ks) catalogs 
with different magnitude limits. They should be used as 
independent catalogs with different magnitude limits and 
areas. This solution can be easily applied to all the cases 
where volume-sensitive statistics have to be applied, as in 
the case of luminosity densities or luminosity functions. 
Analyses that are sensitive to the position information (as 
clustering) instead require the use of the "magnitude limit 
maps" . In the publicly released catalog we provide the cor- 
responding magnitude limit both in z as in Ks with each 
object. 

3. Object detection 

3.1. ACS detection 

The z band of AGS GOODS frames was used as a detec- 
tion image to build a catalog for the GOODS-South field. 
The advantages are twofold: first, the morphological de- 
tails and the resolution of the space-based images cannot 
be obtained with currently available ground-based instru- 
ments in the optical. Second, the longest wavelength avail- 
able in the AGS instrument makes it possible to detect 
high redshift objects (i-dropout). Although fainter sources 
in the i band exist that are not or barely detected in the z 
band images, the significance of their detection is less than 
90 % complete, so they cannot be used as a statistically 
complete sample and are not considered here. 

The optical catalog of GOODS-South field has been 
produced running SExtractor IIBertin fc Arnouts 199611 
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Turning again to the global SExtractor configurations, 
we adopted the parameters listed in Table 01 The adopted 
value for DETECT_THRESH corresponds to a 3.5cr de- 
tection over an area equal to 13 pixels. The high value 
for the cleaning parameter is due to the need to have 
connected pixels in the SEGMENTATION image: with a 
lower cleaning parameter it happens that a segmentation 
with the same identification number can be disconnected 
(not contiguous pixels). 

This parameter set was derived after detailed simula- 
tions and inspections of the background-subtracted and 
object-subtracted images, in order to maximise the com- 
pleteness and to minimise the number of spurious objects 
or artifacts during the detection process. We note in par- 
ticular that wide boxes for the background estimations 
are required to prevent large bright sources to distort the 
background map. 

We compared our catalog with version rl.l of the ACS 
multi-band source catalog released by the GOODS Team 
IjGiavalisco et al. 2004|l . Excluding objects that are fainter 
than the detection limit, there are few galaxies with a dif- 
ference of 1 magnitude or more between our z band esti- 
mation and the one provided by Giavalisco et al. (2004). 
These are systematically blended objects where we used 
the sub-frames with adaptive parameters to optimise the 
detection. Considering only isolated objects, the compar- 
ison gives < \z — Zri.i\ >= 0.01 with a = 0.12. This result 
is not surprising, since both catalogs use the AUTO mag- 
nitudes provided by SExtractor as an estimate of the total 
magnitude. We made no attempt to correct for the known 
biases of the AUTO (or "Kron" ) magnitudes. 

3.2. Simulating the completeness of the detection 

We used a typical value for the r.m.s. in the z band of 
0.0014 (corresponding to a limiting magnitude of z = 
28.18 AB at 1 sigma in 1 sq. arcsec.) to simulate the com- 
pleteness of our detection criterion as a function of the 
magnitude of the sources and their size (half light radius). 
Note that z = 28.18 as a magnitude limit corresponds to 
the deepest exposures for the z band; for the less exposed 
areas, the completeness magnitude scales down to brighter 
values. 

Figure El shows the 90% completeness level resulting 
from such simulations in the z band, for both elliptical and 
spiral galaxies of different half-light radii and bulge/disk 
ratios. Comparing these selection functions with the ob- 
served distribution of magnitude and size for real objects 
in the GOODS field, we derived the exact value of the to- 
tal magnitude z at which our catalog is complete. Table El 
provides the total magnitudes to build a complete catalog 
of sources in the GOODS area. 

The same exercise can be done in the Ks band im- 
ages, where there are 6 main areas, with a broader range 
of magnitude limits. A typical value of 0.06 was used as 
the r.m.s. value for the simulations in the Ks band, corre- 
sponding to a limiting magnitude of 27.0 in a 1 sq. arcsec 



area at 1 sigma. This case corresponds to an exposure time 
of 7 hours, and defines the complete zone at a magnitude 
of Ks = 23.8. 

4. Colour estimate 

As described in Sect. 2, the GOODS data span a 
wide range of image qualities, both in terms of resolu- 
tion/sampling and in terms of limiting depth. In this case, 
it is difficult to design the optimal strategy. On the one 
hand, one would like to make full use of the spatial and 
morphological information contained in the highest qual- 
ity images, limiting the loss of information occurring when 
lower resolution images are included as much as possible. 
This requirement would call for restricting the area over 
which photometry is done. On the other hand, an unbiased 
estimate of the colours is essential for most of the scientific 
application of the GOODS dataset (e.g. photometric red- 
shifts or analysis of the multiwavelength spectral distribu- 
tion). Since most of the detected objects are galaxies that 
may exhibit colour gradients due to a change in the mor- 
phological properties across the wavelengths, one would 
like to ensure that a properly wide area is taken for each 
object. Finally, an optimal S/N is required to improve the 
accuracy of photometric redshifts. 

To match these discordant requirements as closely as 
possible, we adopted two different techniques to accurately 
measure colours, which we adopt in the ACS and in the 
ground-based images, respectively. These methods are de- 
scribed in this section. 

4.1. Colours in the ACS images 

Colours in the ACS images were computed by using the 
isophotal magnitudes as obtained by SExtractor in "dual 
image" mode. This procedure ensures that the photome- 
try is computed at the different wavelengths in the same 
physical region, defined by the isophotal area in the "de- 
tection" image (here, the z band). In the high resolu- 
tion ACS images, the isophotal area follows the appar- 
ent size of the objects and is relatively insensitive to 
the effects of nearby contaminants. For these reasons, it 
has also been used by a number of works dealing with 
similar data set IjLabbe et al. 20031 IVanzella et al. 2001L 
ICimatti et al. 2002|l . 

However, we also explored alternative solutions, all 
based on the "dual image" mode of SExtractor. A fixed 
aperture that is commonly used in the ground-based pho- 
tometry of faint objects is inapplicable for the large dy- 
namical range of the GOODS images: the typical isophotal 
radii range from 2 to 20 pixels, but can be as large as 150 
pixels. 

"Kron" magnitudes ("MAG_AUTO" in SExtractor 
slang) are typically estimated on areas that are larger than 
isophotal ones, and might be less sensitive to colour gra- 
dients. We verified on bright, as well as on faint, isolated 
objects that they provide colours that are (on average) 
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identical to isophotal ones. However, they may be con- 
taminated by neighbouring objects and typically have a 
lower S/N than isophotal, since they include more pixels 
of low S/N. 

We also attempted to use an "optimal" circular aper- 
ture, defined as the aperture which maximises the S/N 
of each object in the detection image. Although promis- 
ing, we verified that it is difficult to provide a robust es- 
timate of such an aperture, since a nearby galaxy may 
easily contaminate the automatic estimate of the S /N of 
faint galaxies. Eventually, we compared the quality of the 
photometric redshifts obtained with different colour esti- 
mators, applied to objects with a spectroscopic redshift, 
finding that isophotal colours provide a slightly better re- 
sult than other choices. On the basis of these results, we 
adopted isophotal colours for all the objects in the sample, 
on all the ACS images. For ground-based images, we de- 
veloped a more complex technique that we discuss below. 



defined by SExtractor. Since the latter typically under- 
estimates the actual object size, the SExtractor isophotal 
area is expanded by an amount that is proportional to the 
object size. 

b) Each object is individually filtered to match the "mea- 
sure" PSF and normalised to unit total fiux: we refer to 
the resulting thumbnails as the "model profiles" of the ob- 
jects. 

c) The intensity of each "model" object is then scaled in 
order to match the intensity of the object in the "mea- 
sure" image. The free parameter for this scaling (Fi) is 
computed with an minimisation over all the pixels of 
the images, and all objects are fitted simultaneously to 
take the effects of blending into account between nearby 
objects. Although the number of free parameters that is 
equal to the number of identified objects is quite large, the 
resulting linear system is very sparse and can be efficiently 
solved using standard numerical techniques. 



4.2. PSF-matched colours in ground-based images 4.2.2. Systematics in the PSF-matching 



To measure colours between ACS and ground-based im- 
ages, we specifically developed and adopted a "PSF- 
matching" code designed to work especially for faint galax- 
ies, which makes it possible to accurately measure colours 
in relatively crowded fields, while making full use of the 
spatial and morphological information contained in the 
highest quality images. 

The technique that we adopted has been introduced 
for the first time by the Stony-Brook group to opti- 
mise the analysis of the J, H, and K images of the 
Hubble Deep Field North. The method is described in 



Fernandez-Soto et al. (1999)] and the catalog obtained 



has been used in several scientific analyses of the HDFN, 
both by the Stony-Brook group fe.g. lLanzetta et al. 19991 
P^hillip ps et al. 2000| ), as well as by several other groups. 
More recently, a conceptually similar method has been 



developed by Papovich et al. (2001) to deal with similar 
problems existing in the first release of GOODS data set. 

4.2.1. The ConvPhot algorithm 

A full description of our code that we name ConvPhot 
is beyond the aim of the present work. Since we make it 
publicly available, we refer the reader to its accompanying 
manual (De Santis et al., in preparation) for a full descrip- 
tion. We recall the basic principle here and highlight the 
possible systematic effects, the strategy that we adopted 
to minimise them, and the validation tests we performed. 

Conceptually, the method is quite straitforward and 
can be better understood by looking at Fig. where we 
plot the case of two objects that are clearly detected in the 
"detection image" but severely blended in the "measure" 
one. The procedure followed by ConvPhot consists of the 
following steps: 

a) Each object is extracted from the "detection image" ^ 
making use of the parameters and of the isophotal area 



As can be appreciated from the example plotted in Fig. 
the main advantage of the method is that it relics on 
the accurate spatial and morphological information con- 
tained in the "detection" image to measure colours in rel- 
atively crowded fields, even in the case where the colours 
of blended objects are markedly different. 

Still, the method relies on a few assumptions that must 
be well understood and taken into account. In particular, 
morphology and positions of the objects should not change 
significantly between the two bandwidths. Also, the depth 
and central bandpass of the "detection" image must ensure 
that most of the objects detected in the "measure" image 
are contained in the catalog. The objects that are deeply 
blended in the measure image should be separated well on 
the detection image. 

In practice, it is unlikely that all these conditions are 
satisfied in real cases. In the case of the match between 
ACS and ground-based Ks images, for instance, very red 
objects may be detected in the Ks band with no counter- 
part in the optical images, and some morphological change 
is expected due to the increasing contribution of the bulge 
in the near-IR bands. Also, in the case where the pixel size 
of the "measure" image (i.e. ISAAC or VIMOS) is much 
larger than in the "detection" one (ACS z in our case), the 
actual limit intrinsic accuracy in image aligning may lead 
to systematic errors. 

To characterise and deal with these possible uncertain- 
ties, we included several options and fine-tuning parame- 
ters in the code to minimise the systematics involved, and 
used extensive simulations to choose them in an optimal 
way. We briefly describe them here, referring to a separate 
paper for a full description. 

Missed flux in the outer area 
The Segmentation image provided by SExtractor is usu- 
ally smaller than the actual object size, since it is limited 
to the threshold used to detect that source. As a result, 
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the object profile and total flux are incorrectly recovered, 
which may result in a systematic bias in the fitted colours. 

To minimise this effect, we expanded the segmentation 
produced by SExtractor by an amount that is proportional 
to the object size, with a minimum area. A minimum value 
for the area of extended segmentation is useful in the ACS 
images of the GOODS-CDFS, where extended but low- 
surface brightness objects are detected only due to the 
bright nucleus, and the isophotal area is limited only to 
the brighter knots. For the GOODS-CDFS z-band image 
we dilate the resulting segmentation such that the objects 
are doubled in linear dimensions, preserving their shape, 
and in the case of very small objects, the area is set to an 
equivalent 1 arcsec (diameter) circular aperture. 

Blended sources in the detection image 
The ConvPhot algorithm is ideally suited to photometry 
of blended objects on the measure image, but requires 
that the same objects in the detection image should be 
well defined and their profiles not be distorted by noise. 
In the practical case, even in the GOODS-ACS images, 
there are galaxies blended in the ACS images or faint ob- 
jects, brighter than the detection limit but still affected by 
noise in their shape/profile. To overcome this problem, we 
introduced an option in ConvPhot to carry out the fitting 
procedure only on the central part of the profile: the fit is 
restricted only to the pixels where the detection image is 
above a given relative threshold. This ensures that the fit 
is carried out only where the signal-to-noise of the model 
(detection image) is high and avoids contamination from 
nearby sources, which cannot be perfectly modelled in the 
detection image. A threshold of 0.5 is used for the GOODS 
dataset, after extensive simulations. 

Alignment errors 
The fitting procedure is extremely sensitive to alignment 
errors. In the simple case the object having a 2-D Gaussian 
shape, it is easy to show that the resulting flux is systemat- 
ically underestimated by a factor / = exp(— 1-^5-), where 
Ar is the offset in pixel of the central position. When 
ground-based images are combined to HST images with 
excellent sampling, this effect is non-negligible. In the case 
of the GOODS ACS data, for instance, the ACS pixel size 
is 0.03", which is often smaller than the residuals of the 
alignment of IR ground-based images. For an alignment 
error of 1 pixel (in the detection image) , a figure that can 
be quite typical or even optimal when combining ground- 
based and HST images, the resulting underestimate may 
be about 3%. 

As a first way out, we included an option in 
ConvPhot to re-centre any object before minimisation. 
In this case, the centre of each object is internally com- 
puted in the model, as well as in the measure image, and 
the measure is re-centered to the model image before min- 
imisation. Since the centre determination may be noisy for 
faint objects, the user can set a limit on the S/N of the 
objects (in each image) for this operation. We use the re- 
centering options only for WFI U band images, only for 



sources with S/N > 15 both in model and measure im- 
ages. 

Variable FWHM or object profile 
Another source of uncertainty may result from a variation 
in the object profile from the detection to the measure 
image, i.e. when the profile in the model image is markedly 
different from the real profile in the measure image, after 
the smoothing with the PSF transformation kernel. This 
can be due to either a physical change of the object profile 
(as due, for instance, to a more prominent bulge in the 
IR) or to an incorrect estimate of the PSF transformation 
kernel. In this case, the resulting flux can therefore be 
either under- or over-estimated, depending on the sign of 
the error in the PSF estimate. An error of 10% in the 
object FWHM will result in a 5% error in the output flux. 

The small systematic effects that have been described 
above, or others resulting from different sources, can be 
efficiently corrected by taking the flux in residual image 
into account. For this purpose, we included an option in 
the code to compute the total residual flux contained in 
the segmentation area of the original frame. This ensures 
the systematics of colour estimations to lower down, as 
shown in the next paragraph and in the ConvPhot paper 
(De Santis et al., in preparation). 



4.2.3. Validation tests 

We performed several validation tests on the ConvPhot 
code, during the debugging phase to estimate the effi- 
ciency in the correction for systematics. The most obvious 
involved the use of simulated images, with a range of lu- 
minosities, PSF, and morphologies, by which we verified 
that the code is computationally correct (De Santis et al., 
in preparation). 

Simulations, however, may not fully reproduce the 
complexity of real objects and data. To obtain a more 
stringent and independent test, we made use of the z band 
FORS image of the CDFS obtained within the K20 survey 
ijCimatti et al. 20T)2| . Here, we used ConvPhot to obtain 
a new estimate of the zacs — zfobs colour in the FORS 
image of the K20. Since the two z filters are fairly simi- 
lar, all objects should have a null colour, barring variable 
objects. In practice a small offset (0.035 mag) between 
the two bands persists still even at bright magnitudes, 
probably due to a different response for the ACS and 
FORS instruments. We limit the comparison to z < 24.75, 
the hmit of FORS z band at S/N = 10. Due to the 
brighter magnitude limit of the FORS image, the error 
in the zacs — zfors colour is dominated by the error 
on the FORS z-band magnitude estimate. The results of 
this comparison are shown in Fig. 13 and show that the 
ConvPhot software is not biased in the colour determina- 
tion. 
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4.3. The final catalogs 

The direct output of a ConvPhot run is the scaling pa- 
rameter Fi for each object. Since the model profile Pi 
for each object is normalised to unit flux, the result- 
ing total magnitude in the measure image is simply 
— 2.51og(i^i) -I- ZPjn, where ZPj^ is the zeropoint of the 
measure image itself. Based on our tests, we concluded 
that this total magnitude is a reliable measure of the ac- 
tual total flux of the objects, somewhat less prone to sys- 
tematic effects than the Kron magnitudes computed by 
SExtractor. However, these total magnitudes can hardly 
be compared with the SExtractor magnitudes of the de- 
tection image, so that reliable colours cannot be obtained 
directly. To this end, we used the total flux Di mea- 
sured by ConvPhot itself in the detection image and 
used it to normahse the object proflle Pi. The resulting 
flux ratio is therefore Flux{measure) / Flux{Detection) = 
Fi/D, X iQ-O-HZP^-zPa)^ jg ^Yie zeropoint of 

the detection image. This flux ratio, or the equivalent mag- 
nitude colour, mmeasure - m detection = "2.5 log(Fi/A) + 

ZPm — ZPd is the final colour that we used in the fol- 
lowing. All colours were finally normalised to the total 
z magnitude to obtain self-consistent magnitudes at all 
wavelengths. 

As stated above, the final colours were obtained with 
ConvPhot. For the J, H , and Ks images, we resampled 
and aligned each individual VLT-ISAAC image to the 
ACS z mosaic that we used for object detection, and used 
the z image, smoothed with the appropriate convolution 
kernel, as the "model" image. In the case of the J/as, L^ag, 
and UviMOS images, we followed the same procedure, 
after cutting the original U images into smaller pieces 
(since ConvPhot can operate only with small size images, 
due to memory limitations), but using the ACS B image, 
smoothed with the appropriate convolution kernel, as a 
"model" to minimise the effects of wavelength-dependent 
morphologies. Because of the geometry of IR images, sev- 
eral objects were detected in more than one IR pointing. 
In these cases, an optimally weighted average was used to 
estimate the final colour. In all cases, we conservatively 
dilated the "Segmentation" image by a factor of 4, with a 
minimum area after the dilate of 800 pixels (corresponding 
to a PSF with 1 arcsec of diameter), and applied a fitting 
threshold of 0.5 to minimise the effects of nearby contami- 
nants. For ISAAC and VIMOS images, we did not use the 
centering option, while for the U35 and t/38 photometry 
we applied a centering with ConvPhot for objects with 
S/N > 15. 

The same technique was applied to the Spitzer images. 
We are aware that the properties of the PSF of the IRAC 
instrument is not fully characterised (FWHM variation 
across the FoV), and the much lower image quality (the 
FWHM ranges from 1.6 to 2.0 arcsec, see Tablenj makes 
the whole analysis more difficult and uncertain. For this 
reason, we regard the colour estimate in the Spitzer bands 
as more uncertain than in the other wavelengths, and we 



plan to address it in more detail in a future paper where 
we will also deal with the Spitzer-selected samples. 

For the Spitzer bands, we used ConvPhot to carry out 
the photometric analysis of the GOODS sources detected 
in the z band. We dilated the "Segmentation" image by 
a factor of 4, with a minimum area after the dilate of 800 
pixels, and applied a fitting threshold of 0.5 to minimise 
the effects of nearby contaminants. For IRAC images, we 
haven't used the centering option. 

Figure [S] shows the colours obtained for all galaxies 
with spectroscopic redshift in our sample (the sources of 
the spectroscopic data adopted are listed in Sect. 6) in 
a few se lected bands. The correspo nding range obtained 
from the Bruzual & Chariot (2003) program is also shown 



for two extreme models, a very young, star-forming model, 
and a maximally old model with a short e-folding star- 
formation timescale. Most of the objects lie within these 
boundaries at the different redshifts, providing a reassur- 
ing quality check on the overall photometry. We verified 
that such an agreement still holds in the other bands. 

5. The Ks-selected catalog 

The catalog produced using the z band of ACS as the 
detection image is not complete in the Ks band magnitude. 
In principle, this could be obtained by using the Ks images 
as detection images, and computing colours in the other 
bands. The catalog obtained, however, would not benefit 
from the higher spatial resolution of the ACS images, and 
the two catalogs would not be homogeneous. To prevent 
this, we followed a different approach of identifying the 
objects missed in the z-selected catalog in the so-called 
"Drop" images produced by ConvPhot in the Ks band, 
i.e. the residuals of the fitting procedure. 

An example is shown in Fig.El which shows that there 
are objects that are not detected in the z band but bright 
in the Ks band. To identify all these missed objects, we 
first carried out a run with ConvPhot with a dilate param- 
eter equal to and using small smoothing kernels in order 
to fit only the central part of the z detected objects and 
not to lose Ks bright galaxies close to z detected objects. 
In this way the Drop images have a minimal area cov- 
ered by "drops" , and the completeness of the Ks selected 
objects is enhanced. 

We ran SExtractor on the Ks "Drop" images, produc- 
ing a complete catalog, whose depth will follow the ex- 
posure maps described above. On the basis of the simu- 
lations described above, it is typically complete down to 
Ks = 23.8 magnitude (AB) at the 90% completeness level. 
As total Ks band magnitude, we used the BEST magni- 
tude of SExtractor. 

The extraction of the colour information for these Ks 
selected objects is not trivial. We calculated the colours 
in an isophotal area defined by the dilated segmentation 
of the Ks band objects. In order to provide an unbiased 
estimate for the colours, we resampled the ACS images to 
the Ks resolution using the convolution kernel adopted by 
ConvPhot. The total magnitude in these bands is com- 
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puted using the Ks as reference and the colour term as 
an additive effect: Bfot = Kstot + {B — Ks)iso, where 
[B — Ks)iso is the colour computed in the "extended" 
isophotal area. For the other bands (U and IRAC) , where 
the FWHM is larger than Ks, we resampled the Ks band 
image to match the measure image PSF. 

At the end we detected 196 objects not present in the 
z selected catalog. The nature of these objects is typically 
twofold: they can be low-surface brightness galaxies in the 
ACS z band images, extended objects that escape from 
detection criterion in z, but with bright total magnitudes 
detected in the Ks band because of a lucky combination 
of pixel-size and seeing effects (151, 77%). The other type 
is constituted by objects that are not detected in the ACS 
image because they are really z-dropout objects, basically 
EROs at z > 1 - 2 (42, 21%). Only three objects are not 
part of these two categories: bright objects in the z and 
Ks bands that are not detected in z because they are on a 
spike of a bright star (2 objects) or in the halo of a bright 
and extended galaxy (1). We recover them from the Ks 
Drop images. 

We merged the z-selected catalog with the Drop Ks se- 
lected objects in order to build a total catalog of objects 
that can be considered complete down to z ~ 26.0 and 
Ks ~ 24.0 simultaneously In the catalog, ID < 20000 
indicates z-detected objects, while ID > 30000 points 
to "drop" Ks selected objects. We also computed, as de- 
scribed in §2.4, the limiting magnitudes (at 1 a and in 1 
sq. arcsec. area) in z and Ks bands for the Ks detected ob- 
jects, in order to have the same quantities of the z-detected 
catalog to ensure a given homogeneity to the total catalog. 
At the moment there is no spectroscopic information for 
the Ks drop detected galaxies, since they are not bright 
enough (in surface brightness limited samples) in the op- 
tical bands or not present in previous ACS-based catalogs 
used for spectroscopic identifications. 

6. The spectroscopic catalog 

The GOODS-CDFS field has been the target of a num- 
ber of spectroscopic surveys that paved this area with 
unprecedented accuracy. The various catalogs surveyed 
this field for different reasons: the COMBO-17 survey 
l|Wolf et al. 2001|l needed spectra to calibrate their pho- 
tometric redshifts over a wide and relatively shallow 
area; the CXO survey ( |Szokoly et al. 2004| ) is a spectro- 
scopic follow-up of X-ray sources in the Megaseconds 
CDFS; and the K20 survey ( |Mignoli et al. 20 05') pro- 
vides the spectra for 327 sources with K < 20{Vega). 
After the public release of the ACS GOODS images, two 
spectroscopic surveys were released: the GOODS VI. 
spectroscopy l|Vanzella et al. 20'05)l and the VVDS survey 
IjLe Fevre et al. 2004|l with the aim of providing a com- 
prehensive census of bright objects in the GOODS field. 
Recently the ESO-GOODS team released a Master'' cat- 
alog of spectra in the GOODS-CDFS area, which sum- 

http://www.eso.org/science/goods/spectroscopy/CDFS_M 



maries all the redshift information of the previously cited 
catalogs, with additional redshifts from sparse surveys. 

The information provided by these surveys is basically 
coordinates and redshifts. The classification (Star, galaxy, 
or AGN) and a quality flag for the spectrum are provided 
by all of the catalogs, but are not present in CXO and 
Master catalogs. The quality flags of the various catalogs 
are not homogeneous, so we tried to define a homogeneous 
classification for all the catalogs. For the K20 catalog, the 
quality of the redshifts is qK20 — 1 for the secure redshifts 
and qK20 = for the uncertain ones. For the GOODS 
spectroscopy, the quality flag goes from A to C, towards 
a decreasing S/N of the spectra. The VVDS uses a re- 
verse scheme when compared to our way of quantifying 
the quality of the spectra, with qvvDS = 4 for the high- 
est quahty to qvvDS = 1 for the lowest one. We define 
four classes of quality flags, from to 3, where qz=0 in- 
dicates the best spectra, with secure identification and 
qz=3 the most noisy spectra. For the other spectroscopic 
catalogs, where the quality flag is not available, a qz=l 
is defined a priori. Only for the K20 and GOODS spec- 
troscopic surveys is there the distinction between early 
or emission line galaxy (or a combined spectrum). For 
the other catalogs, the classification is not homogeneous 
and is not subdivided in finer classifications such as early, 
emission line, or composite spectra. In order to have a sim- 
ple classification scheme, we divided all the objects with 
spectroscopic information in these classes: STAR, AGN, 
GALAXY, EARLY, EMISSION, or COMPOSITE. The 
classification EARLY, EMISSION, or COMPOSITE only 
comes from the K20 and GOODS spectroscopic surveys. 

Since in several cases there are multiple identifications 
for the same object, we cross-correlated the catalogs to 
check for inconsistencies; if there are two different iden- 
tifications for the same object, we take the identification 
with best quality flag. In the GOODS ACS area there 
are currently 1068 known spectroscopic redshifts: Fig. [TUl 
shows the redshift distribution of these objects, in which 
three distinct groups/z-peak are visible, corresponding to 
different groups or sheets of galaxies. There are in total 
928 galaxies, 72 stars, 68 AGNs, and QSOs that are al- 
ready known. In the galaxy sub-sample, only 668 have a 
secure redshift, while the remaining are not secure. For 
428 galaxies the finer classification (early-emission line) is 
available, where "only" 94 are classified as early and 67 as 
composite spectrum. 

The identification of the spectroscopic catalog with ob- 
jects in the ACS images is not trivial, due to the finest res- 
olution of HST compared to ground-based images, where 
close-by objects are merged into a single blob, or because 
there is a coordinate mismatch, especially for the K20 and 
CXO surveys, based on imaging material precedent to the 
ACS one. To overcome this problem/issue, we carried out 
a cross correlation in right ascension and declination be- 
tween the photometric data and the spectroscopic cata- 
log with a relatively large matching radius (1.2 arcsec) 
and divided the objects into two groups: isolated objects, 
t&itht^ unique identification and blended or ambiguous ob- 
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jects, where more than one photometric ID is associated 
with a spectroscopic data. In the last case the associa- 
tion is checked by eye and decided also using the magni- 
tude/colour information. 

Overall, less than 10 percent of the GOODS catalog 
has any spectroscopic information (1068 out of 14847). For 
the remaining objects, having a high photometric quality 
and a wide wavelength sampling of the SEDs, we chose 
the approach of the photometric (§ 7) and neural network 
redshifts (Vanzella et al., in preparation). 

6. 1 . Star-Galaxy separa tion 

We distinguished galaxies from stars and AGNs using 
morphological and photometric information, when spec- 
troscopic data were not available. First, point-like sources 
were selected using the star/galaxy separation flag (s/g) 
provided by SExtractor in the z band. We tuned the se- 
lection on already known spectroscopic stars: objects with 
magnitude z < 20 and s/g > 0.80, 20 < z < 23 and 
s/g > 0.96 or 23 < z < 24 and s/g > 0.98 were con- 
sidered possible AGNs or stars. We used photometric in- 
formation to check this criterion, using in particular the 
"BzK" colour criteria of [Daddi et al. (2004)| Figure [TTl 



shows the positions of AGNs, galaxies, and stars in the B-z 
versus z-Ks two-colour plot. Typically, the point-like ob- 
jects selected with the SExtractor parameter were found 
to lie on the lower side of the diagram, where stars are ex- 
pected to be found. Other objects located in (or nearby) 
this region, but not selected morphologically, were visu- 
ally inspected. Most of them turned out to be extended, 
low-surface brightness galaxies, most likely at very low 
redshift (and hence with similar colour to Galactic stars). 
The remaining unresolved objects were classified as stars. 

7. Photometric redshifts 

On the multicolour catalog described above we applied our 
photometric redshift technique, which was described in 
detail inl Giallongo et al. (1998), and .Fontana et al. (2000) 



(FOO hereafter) which proved to be extremely successful 
in th e HDFN/S p i^ontana et al. 2003) a nd in th e K20 data 
sets (|Gimatti et al. 2002 Fontana _et al. 2004|l . 

In essence, a spectral library of galaxies at arbitrary 
redshifts is computed and a ^^-minimisation procedure is 
applied to find the best -fitting spectral template to the 
observed colours. For each template t at any redshift z we 
first minimise 



E 



observed,! 



s-Ff 



template,! 



(2) 



with respect to the scaling factor s where Fobserved,i is 
the flux observed in a given filter i, (jj is its uncertainty, 
Ftempiate,i is the flux of the template in the same filter, 
and the sum is over the filters used. We then identify 
the best-fitting solution with the lowest Xt z ■ We refer 
to the discussion in FOO for its detailed handling of the 



non-dctcction at faint fluxes. The scaling factor s is there- 
fore applied to the input parameters of the best-fitting 
spectrum to compute all the rest-frame quantities, such 
as absolute magnitudes or stellar masses. 

We tested different choices for the spectral libraries, 
using both empirical templates (based on the Coleman 
et al. data set), as well as on synthetic models, namely 
the iBruzual fc Chariot (2^0 3) and the PEGASE 2.0 
( |Fioc fc Rocca-Volmerang~1997| ). We found that the use 
of synthetic models results in a better agreement with 
the spectroscopic sample available and that, in particu- 
lar, the choice of the PEGASE 2.0 models minimises both 
the r.m.s. and the number of outliers with respect to the 
BC03 library. In this paper, we briefly describe the results 
obtained with the PEGASE 2.0 spectral library to provide 
a quantitative description of the quality of our results, and 
we defer to a forthcoming paper a more accurate compari- 
son of the performances of these libraries, which compares 
these recipes also with the alternative approach of neural 
networks (Vanzell a et al. 2004|l . 

An apphcation of the PEGASE 2.0 code 
to photometric redshifts has been presented by 
Le Borgne fc Rocca-Volmerange (2002)| on the HDFN 



dataset. A major advantage of the PEGASE 2.0 code 
is that it allows to follow the metallicity evolution 
explicitly, also including a self-consistent treatment 
of dust extinction and nebular emission. We adopt 
the same parameter grid as in Cimatti et al. 20021 
and IFontana et al. 20031 which is very close to the 



Le Borgne & Rocca-Volmerange (2002) prescriptions 



In detail, we parameterize the star-formation history 
by two e — folding star- formation time-scales, one (xg) 
describing the time-scale for the gas infall on the galaxy 
and the other (r* ) the efficiency of gas to star conversion. 
By tuning the two time-scales it is possible to reproduce 
a wide range of spectral templates from early types (by 
using small values of Tg and r*) to late. For the earliest 
spectral type, a stellar wind is also assumed to block 
any star- formation activity at an age twind- Dust content 
is followed over the galaxy history as a function of the 
on-going star-formation rate, and an appropriate average 
over possible orientations is computed. The range of 
values for Tg, r*, iujmd and the extinction geometry are 



listed in Table H for a |Rana fc Basu (1992)| IMF. For 
all these models we have adopted a primordial initial 
metallicity. Another differences in our approach compared 
to Le Borgne & Rocca-Volmerange (2002) are that we do 



not apply any constraint on the galaxy ages (apart from 
those set by the Hubble time at each z) and that we have 
allowed for a finer time sampling in the galaxy ages by 
halving the time step with respect to the default value of 
the PEGASE 2.0 program. 

We included the contribution of the nebular emission 
in the synthetic library, both in the continuum and in the 
lines, as allowed by PEGASE 2.0. At the Lya frequency, 
however, inclusion of the corresponding emission line in 
star-forming galaxies would introduce a systematic bias, 
since the observed samples of Lyman break galaxies at 
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Table 4. Parameters used for the PEGASE 2.0 library. 
Other models were added to reproduce the colours of red 
galaxies at high redshift, as described in the text. 



Tg{Myrs) 






extinctiontype 


100 


100 


3000 


spheroid 




100 


500 


5000 


spheroid 




500 


1500 




disk, incl 


— averaged 


1000 


2500 




disk, incl 


— averaged 


1000 


5000 




disk, incl 


— averaged 


2000 


10000 




disk, incl 


— averaged 


2000 


20000 




disk, incl 


— averaged 


5000 


20000 




disk, incl 


— averaged 



z ~ 3 have a wide distribution of observed Lya equivalent 
widths, ranging from absorption to emission systems. For 
this reason, we did not apply the Lya emission line as a 
simple tradeoff, and we also removed the strong absorp- 
tion feature resulting from the stellar features from the 
synthetic spectrum. 

The major difference with the 

approach, 



Le Borgne & Rocca-Volmerange (2002) approach, how- 
ever, is due to evidence that the spectral library described 
above is not able to reproduce the reddest objects (EROs 
or similar) that are typically detected at z > 1 and that 
are either passively evolving galaxies or star-forming 
dusty objects (Cimatti et al 2002) and that we detect in 
our sample. For this reason, we supplemented the above 
models with a few templates designed to mimic these 
objects. 

Passively evolving galaxies were extracted from simple 
truncated models, where a starburst with constant star- 
formation rate is halted after 0.1, 0.3, 0.6, or 1 Gyr. These 
models were used only after the truncation age, so they 
rapidly become redder than the other models. These mod- 
els are able to reproduce redder colours than the stan- 
dard models of PEGASE 2.0, particularly in the near and 
medium infrared part of the SED, thus fitting the ISAAC 
and IRAC observations better. 

Star-forming, dusty objects are instead taken from 
models with a constant star-formation rate, following a 
Calzetti extinction curve with 0.5 < E{B - V) < 1.1. At 
each redshift, galaxies are allowed any age that is compat- 
ible with the Hubble time at that redshift. 

We also added to all models the Lyman series absorp- 
tion produced by the intergalactic medium IjMadau 1995|) , 
but we introduce here a simpler parameterization that 
uses the DA and DB values computed from a compi- 
lation of data at z < 4 and extrapolated at higher z 
with a simple interpolation to the SDSS data. In par- 
ticular we use DA = [(1. + z)2-66838 ^ iq -2.ir28i] and 
DB = £1^-0.3576+1 ^ ;^o-oo36i6 at z < 4, and DA = 
0.187z- 0.247 and DB = 1.178z- 0.133 at 4.3 < z < 6.3, 
with DA = DB = 1 at greater z. 

The most important difference with respect to previ- 
ous papers arises from the inclusion of the Spitzer bands. 



which (at low and intermediate redshifts) extend the ob- 
servation to the spectral regions longer than Sfim rest 
frame. The galaxy emission is not dominated by the inte- 
grated stellar population at A > 5.5/iTO (,Dale et al. 20051 
ILu et al. 2003|l but by the different flavours of dust emis- 
sion, at least in the case of star-forming galaxies that are 
not included in the spectral libraries considered here. As 
a simple way out, we modified our code: any photometric 
band that is contaminated by the dust reprocessed emis- 
sion above 5.5/iTO in the rest-frame is ignored in the fit. 

We tested our recipe against the large subsamplc of 
galaxies with secure spectroscopic redshift, which wc com- 
piled as described in the previous section: to this end, we 
considered 668 galaxies with the spectroscopic flag equal 
to or 1. The results are shown in Fig ll2l In its upper 
panel, we plot the Zspec—Zptot relation to show that we find 
an excellent agreement between photometric and spec- 
troscopic redshifts over the fully accessible redshift range 
< z < 6, with a very limited number of catastrophic er- 
rors. More quantitatively, we plot the distribution of the 
absolute scatter Az = {zgpec — Zphot) in the inset: the cen- 
tral part of the distribution is perfectly represented by a 
Gaussian distribution (smooth red curve) with a standard 
deviation a — 0.06 and a small number of outliers. Because 
of these outliers, the distribution is not Gaussian and the 
usual r.m.s. is not a good indicator of the width of the dis- 
tribution. If we instead adopt the average absolute scatter, 
we find < |Az| >= 0.08 and < |Az/(l -I- z)| >= 0.045: 
these values are among the lowest ever obtained with 
the photometric redshift technique in the redshift inter- 
val < z < 6. Similarly, the lower panel of Figll2l 



shows the relative scatter (zgpec — Zphot)/{^ 



•^spec 



re- 



stricted to the z < 2 range and discarding the most dis- 
crepant objects, which has a nearly homogeneous r.m.s. 
of 0.03. As far as we know, this is the highest precision 
ever obtained on faint samples spanning such a wide red- 
shift range, although it is still lower than the average 
SDSS value. Several elements concur to obtain this im- 
provement with respect to similar samples (e.g. Fontana 
et al 2003, MoE asher et al. 20041 IVanzella et al. 2004|l : 
the larger number of bands, the sophisticated technique 
adopted for photometry, and the use of PEGASE 2.0 mod- 
els, which provided slightly better results than other syn- 
thetic models or observed templates. We plan to discuss 
these effects better in a future paper. 

The resulting redshift distributions are shown in 
FigEl both for the 9862 galaxies of the z-selected sam- 
ple and for the 2931 galaxies of the Ks-selected one. They 
were drawn by adopting the photometric redshift for all 
galaxies except those with a secure spectroscopic one. For 
the reasons described in Sect. 4, the two samples do not 
have a unique magnitude limit, but for most of the sample 
the typical magnitude limits are Ks < 23.8 and z < 26. 
The two distributions show both a marked peak at z ~ 0.7, 
which is a well known feature of this field, and a marked 
decrease at z > 1.5. 
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8. Summary 

In this paper, we have described in detail the proce- 
dures adopted to obtain a multiwavelength catalog of the 
large and deep areas in the GOODS Southern Field cov- 
ered by the deep near-IR observations obtained with the 
ESO VLT. The catalog, named GOODS-MUSIC (MUlti- 
wavelength Southern Infra-red Catalog) , is made publicly 
available and we plan to use it in several future scientific 
analyses. The main features of this work are the following: 

• The catalog is based entirely on public data: we in- 
cluded the FA35W, i^606W^, F775W and F850LP ACS 
images in our analysis, along with the JHKs VLT data, 
the Spitzer data provided by the IRAC instrument (3.6, 
4.5, 5.8, and 8.0 /im), and publicly available U-band 
data from the 2.2ESO (C/35 and f/gg) and VLT-VIMOS 

{UviMOs)- 

• The complications arising from the inhomogeneous 
coverage and depth of the available images were addressed 
with great care. The first is that a unique magnitude 
limit cannot be defined in a given band. For this purpose, 
we included a careful estimate of the magnitude limits 
across the detection images (z and Ks) to make the ex- 
traction of statistically well-defined samples possible. In 
practice, both the z and the iiTs-selected catalogs were 
split into several sub-catalogs with a well-defined magni- 
tude limit, using the bins listed in Tab[21 These catalogs 
should be used as independent catalogs in all those cases 
where volume-sensitive statistics have to be applied, as in 
the case of luminosity densities or luminosity functions. 
Analyses that are sensitive to the position information 
(such as clustering) instead requires the full use of the 
positional information. 

• The object detection was first done in the ACS z (i.e. 
F850LP) image, using a customised version of SExtractor, 
which we developed to cope with the large dynamical and 
morphological range of the ACS mosaics. In addition to 
this, we identified the objects that are detected in the Ks 
images that escaped detection in the z one. This double- 
pass procedure enabled us to obtain a unique catalog from 
which either a ^-selected or a ifs-selected subsample can 
be easily extracted. 

• A detailed set of simulations was executed to esti- 
mate the completeness limits of the samples. Although the 
effects of inhomogeneous coverage and depth mentioned 
above prevent a unique threshold from being adopted, we 
find that the typical completeness is at the level of z ~ 26 
or Ks ~ 23.8. 

• The galaxy colours in the ACS images were esti- 
mated straightforwardly using the isophotal magnitudes 
provided by SExtractor. Conversely, the colour estimate 
on the IR or UV images, which have much poorer quality, 
was done using a specific "PSF-matching" software that 
we developed to this end and that we named ConvPhot. 
It is described at length in Sect. 4. As a result, all the 
14847 galaxies (plus stars and AGNs) in the sample have 
a 1/35,1/38, B,V,i, z, J, Ks, 3.6, 4.5, 5.8, 8.0 coverage, 9651 
also have UyiMos and 8441 have H. 



• We cross-correlated our catalog with all the spectro- 
scopic surveys available to date, assigning a spectroscopic 
redshift to more than 1000 sources in our catalog. 

• The final catalog is made up of 14847 objects, with 
at least 72 known stars, 68 AGNs, and 928 galaxies with 
spectroscopic redshift (668 galaxies with reliable redshift 
determination) . 

• We applied our photometric redshift code to the 
14 bands catalog. We applied a standard technique, 
choosing a set of synthetic templates drawn from the 
PEGASE2.0 synthesis model. The comparison with the 
spectroscopic sample shows that the quality of the result- 
ing photometric redshifts is excellent, with an r.m.s. scat- 
ter of only 0.06 for the redshift interval < z < 6. 

• The fuU multicolour GOODS-MUSIC catalog, in- 
cluding the redshift information (both spectroscopic and 
photometric) is made publicly available together with the 
software specifically designed for this purpose at the site 
http://lbc.oa-roma.inaf.it/goods or at CDS^. An example 
of the public data can be found in Tables El and H 
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ID 


RA 
J2000 


DEC 
J2000 




class'' 


catalog'^ 




zphot 


POS" 


star* 


AGN^^ 


zlim'' 


Kslim' 


S/Gj 


10015 


53.109565 


-27.788200 


0.995 


emission 


K20 





0.995 


1 








28.044 


26.982 


0.030 


10016 


53.166175 


-27.787519 


1.097 


galaxy 


C0MB017 





1.097 


1 








28.026 


27.008 


0.000 


10017 


53.124374 


-27.788923 


-1.00 


unknown 


unknown 


99 


0.300 


1 








28.036 


26.957 


0.010 


10018 


53.056438 


-27.788972 


-1.00 


unknown 


unknown 


99 


3.500 


1 








28.046 


26.327 


0.010 


10019 


53.046402 


-27.789001 


-1.00 


unknown 


unknown 


99 


0.280 


1 








28.051 


26.136 


0.740 


10020 


53.065872 


-27.787111 


0.738 


early 


K20 





0.738 


1 








28.036 


26.326 


0.030 



[^) spectroscopic rcdshift (-1.0=not available). 
[^) spectroscopic class (see text). 
['^) reference spectroscopic catalog. 

['^) quality of spectroscopic redshift (0= very good, l=good, 2=uncertain, 3=bad quality, 99=not available). 
[") position flag (l=inside GOODS-MUSIC area, 0=outside). 

'*) star flag (Improbable star, 0=no star), on the basis of spectroscopy, morphology, and BzK colours (see text). 

'^) AGN flag, based only on spectroscopy (l=probable AGN, 0=no AGN). A galaxy should have star flag=0 and AGN flag=0. 

[^) magnitude limit in the z band in 1 sq. arcsec. and at 1 a. 

J') magnitude limit in the Ks band in 1 sq. arcsec. and at 1 a. 

[') Star/Galaxy index of SExtractor. 



Table 6. Extract of the GOODS-MUSIC web catalog: optical photometry. 



ID 


Ua5 


U:,8 


UviM 


B 


V 


i 


z 


(735 


t/38 


UviM 


B 


V 


i 


z 




mag 


mag 


mag 


mag 


mag 


mag 


mag 


err 


err 


err 


err 


err 


err 


err 


10015 


23.108 


23.073 


99.0 


23.130 


22.851 


22.204 


21.905 


0.017 


0.044 


99.0 


0.009 


0.006 


0.007 


0.006 


10016 


23.085 


22.954 


99.0 


23.029 


22.618 


21.938 


21.425 


0.031 


0.070 


99.0 


0.017 


0.010 


0.011 


0.008 


10017 


26.356 


-25.417 


99.0 


28.096 


26.397 


25.339 


24.992 


0.965 


1.085 


99.0 


0.207 


0.084 


0.084 


0.067 


10018 


26.806 


24.918 


99.0 


27.257 


25.828 


25.444 


25.523 


0.820 


0.477 


99.0 


0.274 


0.046 


0.062 


0.086 


10019 


27.843 


-26.546 


99.0 


29.091 


26.807 


26.547 


26.550 


1.014 


1.085 


99.0 


0.630 


0.076 


0.107 


0.133 


10020 


24.893 


24.945 


99.0 


24.254 


22.429 


20.959 


20.473 


0.135 


0.386 


99.0 


0.044 


0.007 


0.003 


0.002 



Magnitudes are in the AB photometric system. Negative magnitudes indicate upper limits, while mag = 99.0 and err = 99.0 
indicate that the measure is not available. 
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Table 7. Extract of the GOODS-MUSIC web catalog: infrared photonietry. 



ID 


J 


H 


Ks 


3.6/i 


4.5/i 


5.8/i 


8.0m 


J 


H 


Ks 


3.6/Lt 


4.5/i 


5.8m 


8.0/i 




mag 


mag 


mag 


mag 


mag 


mag 


mag 


err 


err 


err 


err 


err 


err 


err 


10015 


21.739 


21.534 


21.468 


21.523 


21.780 


22.162 


22.671 


0.017 


0.024 


0.013 


0.009 


0.019 


0.095 


0.150 


10016 


20.971 


20.662 


20.235 


20.142 


20.328 


20.793 


20.835 


0.011 


0.013 


0.010 


0.009 


0.010 


0.032 


0.033 


10017 


24.696 


24.265 


24.018 


25.739 


-25.750 


-24.278 


-24.275 


0.237 


0.308 


0.132 


0.509 


1.085 


1.085 


1.085 


10018 


25.579 


-26.234 


25.129 


25.086 


25.071 


23.713 


-23.970 


0.429 


1.085 


0.527 


0.412 


0.691 


0.816 


1.085 


10019 


-26.806 


26.237 


25.769 


24.371 


25.263 


24.543 


23.041 


1.085 


0.905 


0.799 


0.155 


0.348 


0.914 


0.210 


10020 


20.000 


19.496 


19.218 


19.267 


19.722 


20.007 


20.813 


0.005 


0.004 


0.004 


0.003 


0.004 


0.017 


0.034 



Magnitudes are in the AB photometric system. Negative magnitudes indicate upper limits, while mag = 99.0 and err = 99.0 
indicate that the measure is not available. 



Steidel, C. C, Pettini, M., Hamilton, D. 1995, AJ, 110, 2519 
Steidel, C. C, Adelbergcr, K. L., Giavalisco, M., Dickinson, 

M., Pettini, M. 1999, ApJ, 519, 1 
Szokoly, G. P., Bergeron, J., Hasinger, G., Lehmann, I., 

Kcwlcy, L., Mainicri, V., Nonirio, M., Rosati, P., Giacconi, 

R., Gilli, R. et al. 2004, ApJS, 155, 271 
Valdes, F. G. 2002, ADAAS, 309; Edited by Ranjan Gupta, 

Harinder P. Singh, Coryn A.L. Bailer-Jones. New Delhi; 

London 

Vandamc, B. 2002, SPIE, 4847, 123 
Vanzella, E. et al. 2001, AJ, 122, 2190 

Vanzella, E., Cristiani, S., Fontana, A., Nonino, M., Arnouts, 
S., Giallongo, E., Grazian, A., Fasano, G., Popesso, P., 
Saracco, P., Zaggia, S. 2004, A&A, 423, 761 

Vanzella, E., Cristiani, S., Dickinson, M., Kuntschncr, H., 
Moustakas, L. A., Nonino, M., Rosati, P., Stern, D., 
Cesarsky, C, Ettori, S., et al. 2005, A&A, 434, 53 

Wolf, C, Mcisenheimer, K., Roser, H. J., Beckwith, S. V. W., 
Chaffee, F. H. Jr., Fried, J., Hippelein, H., Huang, J. S., 
Kummel, M., von Kuhlmann, B., Maier, C, Phleps, S., 
Rix, H. W., Thommes, E., Thompson, D. 2001, A&:A365, 
681 



Appendix A: True RMS estimation with 
correlation matrix 

Wc provide detailed calculations here concerning the cor- 
relation of the noise and the theoretical interpretation. 

Suppose we have a l-dimensional distribution Xi, with 
N random values. This distribution has a mean m and 
an r.m.s. ax- We define the correlation between adjacent 
pixels of distance equal to j as: 



1 ^ 



(A.1) 



The term Co is the variance of the distribution cr^ . If the 
pixels are all randomly distributed and uncorrelated, the 
general term Cj behaves as a Delta function: 



= 

7^0 



(A.2) 



The terms Cj should be nearly symmetric, if the number 
N is large, C-„ C+„. 

Considering the case of a correlated distribution y built 
from this uncorrelated distribution, we have: 



E+k 
h=-k WhXi+h 



E+k 
h=-k Wh 



(A.3) 



If the kernel w conserves the flux, then ^h=-k'^h = 1 

and yi = J2h=-k'^hXi+h- 

We do not demonstrate (trivial) that the two distribu- 
tion Xi and yi have the same mean = TUy. We consider 
here that the mean is null, m = 0. 

We now want to also evaluate the correlation coeffi- 
cient Cj for the new distribution yf. 



1 ^ 

i=l 

which becomes 



(A.4) 



1 ^ 



N 



+k 

^ WhXi+hVi+j 

h=-k 



i=l 

Now we can define 

N 

N 



+k ^ N 



'N 

h=-k i=l 



1 ^ 

i=l 

which yields 

AT 



Ui 



hi 



or 



-y 



Xi+h 



+k 

^ WlXi+j+l 
l=-k 



+k 
l=-k 



1 ^ 

-T 

N ^ 



^i-\-h'^i-\-j-\-l 



i=l 



j+l-h ■ 



(A.6) 

(A.7) 
(A.8) 

(A.9) 



l=-k 
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Table A.l. Correlation matrix of noise in the ACS z band 



Column 


-1 





+1 


+1 


0.06411 


0.28940 


0.06131 





0.28436 


1.00000 


0.28436 


-1 


0.06131 


0.28940 


0.06411 



Now we can sum over the terms CJ: 



+m 



+fc 



E E '^I'^hi-h 



j — — m \_h— — k l — — k 



+m 



-\-k -\-k +771 

E = E -'^ E E ^^^-^ 

h— — k l — — k j — ~7n 



(A.IO) 



(All) 



'l2j=~rn ^j+l-h 



+m 



Cj+i-h — '^l: ill that case we have 

+k +k +k +k 



E = E ^'^ E = CT^ 5] u;, 5] .(A.12) 



h— — k l— — k 



h— — k l — — k 



Using the fact that the kernel conserves the flux, 



J2ht-k'^h = and Y^'^^^^wi ^ 1, it finally yields 



+ 771 

E 



E 

j=~m 



(A.13) 



Hence, we have demonstrated that the sum of the corre- 
lation coefhcient Cj gives the total variance of a distri- 
bution, whether correlated or not. Indeed, since Cq — ay 

and Cq < J2^=-m '^^''^ conclude that the variance 

of a correlated distribution is less than the variance of the 
original distribution, or true variance, and it is an under- 
estimated value for the r.m.s. of the data. 

These results are derived for a discrete distribution, 
but they can be generalised to the continuum case with the 
help of integrals and Fourier Transform. All these results 
derived from a 1-dimensional distribution can be easily 
generalised for a 2-dimensional distribution, which in our 
case is a reduced image, where the noise between adjacent 
pixels is correlated, due to the algorithms used to sum the 
dithered exposures and the seeing. 

In the z band, the resulting kernel (Tab. lA.ip is posi- 
tive in the central 3x3 pixels, and fades rapidly to zero at 
large distances. If one only considers the central 3x3 pixels 
and normalises the peak value to 1, it results a correlation 
matrix as in Tab. lA.ll 

The sum of all the terms gives the ratio between the 
true and the measured variance (c^/c^) of the image. 




Z BAND ACS 




MAGLIM (AB) 

Fig. 3. The magnitude limits (at 1 sigma and in 1 sq. arc- 
sec) for the z band in the GOODS area. Dark areas cor- 
respond to deeper exposures, while white zones are much 
shallower ones, as also shown in the histogram of the mag- 
nitude limits. The geometry of the dark and white areas 
depends on the dithering strategy adopted. 
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Fig. 4. The magnitude limits (at 1 sigma and in 1 sq. 
arcsec) for the Ks band in the GOODS area. Dark ar- 
eas correspond to deeper exposures, while white zones are 
much shallower ones, as also shown in the histogram of 
the magnitude limits. 
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Fig. 5. The half-light radius of galaxies in the GOODS- 
CDFS region as a function of the observed z and Ks 
magnitudes. The solid line indicates the completeness at 
90% value for elliptical galaxies, while the short-dashed 
line shows the same for spiral galaxies. We simulated 100 
galaxies for each bin of magnitude, half-light radius, and 
bulge/disk ratio: the magnitude range is from 20 to 30 in 
steps of 0.5, the half-light radius range from 0.1 to 1.5 in 
steps of 0.1 arcsec, and the B/D from 0.3 to 0.9 in steps of 
0.1 for ellipticals or from 0.02 to 0.08 in steps of 0.01 for 
spiral galaxies. The completeness as a function of magni- 
tudes and half-light radii is the fraction of galaxies recov- 
ered by our detection procedure, averaged over the whole 
bulge/disk ratio interval. Long-dashed lines represent the 
fainter envelope of the observed galaxies. The complete- 



A. Grazian et al.; GOODS- 



a; 



b) 




Fig. 6. A schematic representation of the ConvPhot algo- 
rithm. 

a). Two objects are clearly detected and separated in the 
high resolution detection image (blue, solid-thin line). The 
same two objects are blended in the low resolution ^^mea- 
sur^^ image (red, solid-thick line) and have quite different 
colours, h). The two objects are isolated in the high reso- 
lution detection image and are individually smoothed to 
the PSF of the measure image to obtain the "model" im- 
ages, c). The intensity of each object is scaled to match 
the global profile of the measure image. The scaling factors 
are found with a global minimisation over the whole 
object pixels. 
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Fig. 7. The zacs — Zfors colour as a function of the mag- 
nitude of the objects detected on the K20 survey (upper 
panel). The histogram of this colour (lower panel) shows 
that the ConvPhot magnitude determination is not biased 
and is compatible with aperture photometry. 
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Fig. 9. A so-called "Drop" image produced by ConvPhot 
in the Ks band. The name derives from the fact that each 
z-band detected object that is fitted by ConvPhot is mul- 
tiplied by zero (drop) in order to leave only those Ks bright 
galaxies not detected in the z band. 



Fig. 8. Observed colours as a function of redshift for 
the objects with spectroscopic redshifts in the GOODS- 
S field. The solid lines are (from top to bottom) the 
colours predicted by the Bruzual & Chariot (2003) code 
for the following cases: (i) a maximally old model, ob- 
tained assuming an exponential history of star-formation 
with timescale 0.1 Gyrs started at z = 10 with metallicity 
Z = 2.5Zq; (ii) a similar model but with an exponential 
timescale of 5 Gyrs and an actively star-forming and non- 
evolving model with a constant star-formation rate and 
age of 1 Gyr. Dotted lines represent the redshift range in 
which the colours are not reliable due to the incomplete 
modelling of dust emission at A > 5.5/im rest frame (see 
Sect. 7). 




Fig. 10. The histogram of spectroscopic redshifts for 
galaxies in the GOODS-CDFS field. Three peaks in red- 
shifts show the presence of large-scale structures like 
groups or sheet /wall of galaxies. 
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n AGN 

J- galaxies 1.4<z<3.5 
• spectroscopic stars 

-■*■- point lilve sources 



4 6 8 

(n-/) 

Fig. 11. The B-z versus z-Ks colours for objects in the 
GOODS south. Small dots are galaxies with redshifts z < 
1.4 or z > 2.5, squares are spectroscopically known AGNs, 
arrows indicate galaxies with 1.4 < z < 2.5 (spectroscopic 
and photometric redshifts), filled points are spectroscopi- 
cally known stars, while stars represent objects with point- 
like morphology according to SExtractor. Solid lines in- 
dicate the BzK criterion to isolate star-forming (upper 
left part of the diagram) and passively evolving galax- 
ies (upper right part of the diagram) at 1.4 < z < 2.5. 
Objects below the dashed line are stars, according to 



Daddi et al. (2004) 
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Fig. 12. Upper panel: the relation between the spec- 
troscopic (x-axis) and the photometric (y-axis) redshift 
on 668 galaxies with accurate spectroscopic redshift. In 
the inset, the distribution of the absolute scatter Az ~ 
{zspec — Zphot) is shown and compared with a Gaussian dis- 
tribution with a standard deviation a = 0.06 (smooth red 
curve). Lower panel: relative scatter {zspec — Zphot)/i^ + 
Zspec), restricted to the z < 2 range and discarding the 
most discrepant objects in the same sample. 
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Fig. 13. Upper panel: redshift distribution of 9862 galax- 
ies in the z-selected sample. The typical magnitude limit 
is ~ 26. Lower panel: redshift distribution of 2931 galax- 
ies in the Ks-sclected sample. The typical magnitude limit 
is Ks ~ 23.8. 



